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ABSTRACT 



A method is proposed to relate simple biological 
response modeling to air quality criteria. We also suggest 
a procedure using the biological response curve in conjunction 
with measured pollutant concentrations for setting emission 
levels to ensure compliance with air quality standards. A 
byproduct of our analysis is an air pollution index which 
avoids the arbitrariness of commonly used indices. 



1. Introduction 

In a series of papers (see Larsen et al. , 1977) 
Larsen has proposed a rational procedure to set air quality 
criteria. His method is based on an extensive study of the 
response of vegetation and mice to exposure from various 
pollutants. Larsen found that the pollutant concentration 
required to produce a specified effect varied as some power 
of exposure time. In other words, the concentration versus 
exposure duration plotted as a straight line on a log-log 
graph. Clearly this type of plot is useful for the derivation 
of air quality criteria (Munn et al. , 1977). Instead of 
concentrations, (which are often chosen independently of each 
other) at distinct averaging times we have an "operating- 
curve of air quality criteria as a continuous function of 
averaging time. While this approach is preferable to commonly 
used criteria setting procedures it is still a static method 
in that it does not account for the continuous response of 
organisms to pollutant doses. This disadvantage becomes clear 
when we consider a situation where say two hours of concentrations 
slightly below the criterion are followed by one hour in which 
the criterion is violated. Obviously, this one hour violation 
has a greater effect than that preceded by hours of zero con- 
centration. Air quality criteria do not distinguish between 
these two situations. 

In this paper we present a simple dynamic model of 
biological response to pollution. We show that this type of 
model, which has been used by Saltzman (1970) to study bio- 
logical transfer functions J can also be utilized to derive air 
quality criteria. 



- 1 - 



2. Model 

The simplest model of biological response can be 
represented by the equation 

In (1) ot and ft are proportionality constants , and E represents 
the effect produced by the ambient pollutant concentration C a . 
It is clear that we assume that the pollutant uptake rate is 
proportional to C and the detoxification rate is proportional 
to E. The meaning of (1) is clarified if we take C & to be the 
concentration of CO and E to be that of carboxyhemoglobin (COHb) . 
There is some research (Roach, 1966) which indicates that the 
response of COHb in human blood to CO exposure is modeled fairly 
well by the solution of (1). Roach also shows that a wide 
variety of pollutants elicit this type of response from critical 
organs in the human body. We feel that the extension of the 
model to vegetation response is not unreasonable. 

For a constant C ,the solution of (1) is 

We assume that a specified degree of pollutant effect 
corresponds to some definite level of E. Then, taking 
E(o) =0, the concentration required to produce a specified 
effect over a given exposure time t is given by 

n C m (3) 

C at " (l-exp(-=< t)) 

In (3), C is the pollutant concentration which produces the 
specified effect when t^"oO . In other words, C m is the con- 
centration to which the organism can be exposed continuously 
with no "adverse" effects. It is noted that C m essentially 
specifies the degree of injury. As £> is arbitrary, we take 
o<- G> so that the allowable effect becomes C m , the long term 
concentration. For small t, Eq. 3 reduces to 

Ok . -t = C (4a) 

at m 
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Thi:> shows that equal doses produce equal effects only at small exposure 
tim;-. . This can be clarified by rewriting (3) as 

C m = D (1 - exp (-pit ) (4a) 

As C is the effect variable, the effect is proportional to the dose 
in 

D = C . t only for small t. In fact, for a constant dose the effect 

at 

decreases as the exposure time increases. This result is compatible 
with experimental findings (Heck and Brandt, 1977). 
3 . A ir Quality Criteria 

It is clear that Eq. (3) can be used to derive air 
quality criteria. In Fig.l we have plotted (3) on a log-log 
graph. The curve, which we refer to as the biological response 
curve BR, is specified by *t (which is the inverse of the 
biological response time ^i \ and the effect C^. In the 
figure, the straight line D^^ is a representation of measured 
maximum concentrations as a function of averaging time. The 
form of D, is based on several studies in urban areas (Larsen, 
1971) which indicate that concentration versus averaging time 
plots as a straight line on a log-log graph. If we desire 
to prevent biological effects above a level characterized 

by c it is clear that measured concentrations given by D, 

m * 

have to be reduced. Between the averaging times t A and t B 
expected maximum concentrations are above the limits set by 
BR. In order to prevent pollutant induced injury we have to 
shift the curve D, so that it lies below BR at all averaging 
times. The line D„ , which is tangential to BR at C, satisfies 
the air quality criteria set by BR. It is clear from the 
figure that the critical averaging time t emt the associated 
concentration is determined by the condition 

Slope of BR at t = Slope of D^ (= slooe of D 2 ) 

We note that emission reduction does not chanqe the slope of 
the measured concentration. It is readily seen that the 
above condition reduces to 



1*o+fl-l}/(\-£ff(-<fi y y) = j> 



(5) 



^-f- 



W<? 



where v - otu and b is the negative of the slope of D^ 
note that (5) | s independent of C . Thus, once we know x 
as a function of p, we can determine the point through v;hich 
D should pass by specifying the effect C 



Let us illustrate the procedure described in the 

previous paragraph by considering carbon monoxide, a gas 

whose effects on the human body are fairly well known. It 

is recalled that CO standards are based on criteria relevant 

to the human body. Calculations by Roach (1966) indicate 

that ~C. for CO is 2.88 hr. Then, assuming that the 1-hour 

b 
(MOE, Canada) standard for CO of 30 ppm is related to a 

permissible level of COHb, we can compute C which turns out 
to be 8.80 ppm. With the biological response curve completely 
specified we can compute t and the associated CO concentration 
for a specific urban area. For a typical value of p we 
choose 0.15 from a study by McGuire and' Noll (1971) . From 
Fig. 2 we find that for p = 0.15, x is 3.06, a value which 
translates into t = 8.81 hr. From (3) we can compute the 
concentration as 9.23 ppm. The calculation indicates that 
for an area with p = 0.15 we can satisfy the criteria set by 
the BR curve for CO by ensuring that the maximum concentration 
for an averaging time of 8.81 hr. is less than or equal to 
9.23 ppm. It is worthwhile to clarify the meaning of this 
criterion. What we are saying is that by meeting the 8.81 hr. 
concentration requirement we do not have to worry about con- 
centrations at other averaging times. This would be true, of 
course, only if the maximum concentrations did not deviate 
from the D curve. As p is a function of time and area for the 
same pollutant (McGuire and Noll, 1971), it is clearly not very 
meaningful to derive criteria and hence standards from the D 
line. The criteria are set by the BR curve. While 
the D line can be used to help us to meet the requirements 
of the BR curve, it cannot specify a set of objective standards, 

The discussion cf the previous paragraphs shows that 
the concept of the biological response curve can be very use- 
ful in the establishing of air quality criteria. However, its 
utilization forces us to replace time varying concentrations 
by averages. Clearly, the response to a varying concentration 
over a soecified interval of time is not equivalent to that 
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to a constant average concentration over the same time period. 
In the next section, we propose a dynamic model of biological 
response in order to improve upon the "static" nature of the 
derivation of criteria we have just discussed. 

4 . Dynamic Model and An Air Pollution Index 

The general solution of Eq. (1) can be written as 

-t 

As C (t) is known, it is a simple matter to compute E(t). 

a 
Then, from the previous discussion it is clear that the state- 
ment of our "dynamic" criterion is that E(t) should not exceed 

C which specifies the level of allowable injury. This concept 
m 

can be expressed in a diagram. 

Criterion: E(t)< c 



C, (t) 




E (t) 



What we are essentially doing is transforming the ambient 
concentration C (t) to E(t) through the parameter d (= V^u) • 
The integral sign indicates that E(t) represents the cumulative 
effect due to C (f ) . 

CI 

The preceding discussion suggests that we can construct 
an air pollution index by normalizing Eq. (6) with C . Then the 
equation describing the Index I = E(t)/C is: 



1(0 -TtoVfrf t-rfft +" ^L K^T^p^^ p^t)^ t 



(7) 



nvt.* p 



We notice that I ranges from to 1 when the c umulative effect 
of pollutants is within allowable limits. An index greater than 
one indicates a problem. 
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An approximate form of (7) can be written as 




c 



tn. 



(8) 



The attractive feature of the index in (8) is that it 
accounts for the cumulative effect of pollutants. Furthermore, 
as (8) mimicks the actual biological response, I is a measure 
of the injury caused by the pollutant in question. Clearly, 
this is a vast improvement over the commonly available air 
pollution indices (Thorn and Ott,19 75) the basis of which are 
admittedly arbitrary. 

Knowledge of C and o(, for commonly found pollutants 
will allow us to compute indices for each of the pollutants. 
For the purpose of reporting to the public, individual index 
values or the maximum value can be used. It is noted that the 
suggested index does not account for synergistic effects. But 
then, as far as we are aware, no other index does so in a 
scientifically acceptable manner. 

5. Determination ofe/ and C 
_ & rn 

As our model of biological response depends on the 
two parameters*^ and C , in theory, we can specify BR by 
experimentally determining two points on the curve. Suggestions 
on the details of an experimental program to determine these 
parameters are outside the scope of this paper. However, it is 
worthwhile to deduce o< and C from secondary air quality 
standards, assuming of course, that they are "based" on points 
on the BR curve. 

One of the graphs in Fig. 3 is a BR curve for vegetation 
based on U.S.-E.P.A. secondary standards (3-Hr., 24-Hr., and 
annual) for S0_. By choosing two standards at a time we can 
fit curves to generate three sets of C andcx^ (- IjnCJ . At the 
present time there is not enough information to choose between 
the BR curves. However, as the 3-Hr. and 2 4-Hr. concentrations 
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are probably better established than the annual concentration, 

3 
the C value of 180 ug/m and T = 19.8 hrs miqht be the best 
m £, 

parameters for the formulation of an air pollution index. 

Figure 4 shows the BR curves for CO generated using the 

1-hr and 8-hr standards set by two air control agencies, 
the Ministry of the Environment, Ontario, and the United 
States E.P.A. We see that while the 1-hr standard set by 
MOE is lower than that set by E.P.A. , the 8-hr MOE con- 
centration is higher than the E.P.A., standard. This 
difference implies that MOE standards are based on a biological 
response time of 1.79 hrs while the E.P.A. implicity assumes 
that the human body takes a longer time of 3.96 hours to 
recover from CO dosage. We also note that the MOE curve 
indicates that no adverse effects result from a continuous 
exposure of 12.8 ppm of CO. On the other hand, the E.P.A. curve 
fixes this long term concentration at a lower value of 7.81 ppm. 
It is necessary to recall that the preceding interpretation 
of the CO standards is based entirely on the hypothesized 
biological response model. 

In the same figure we have drawn straight lines to 
indicate the consequences of a CO data line which passes 
through the points specified by the standards. We see that 
as the straight lines lie above the BR curves between the 
standard points, meeting the 1-hr and 8-hr standards does not 
ensure that the criteria set by the BR curve' are met. In view 
of this, we feel that the best way of setting CO standards is 
by using the criteria derivation procedure described in section 
3. In fact, for pollutants such as CO which are associated with 
relatively short biological response times, it would be 
necessary to use a dynamic model to stay below the specified 
injury level. 



6 . SUMMARY 



A simple dynamic model of biological response to 
pollutant dosage can be used to derive a plausible form of the 
biological . response curve. -We have shown how the BR curve 
:carTbe used to. establish air- quality criteria. r. We feel- that 
the formalism of the^derivation :is important by itself. . 
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The response model can also form the basis of an 
air pollution index which avoids the arbitrariness of currently 
available indices. The index is derived from sound physical 
principles and is thus preferable to other indices even if 
the parameters C< and C have to be chosen in a less than 
satisfactory manner. 

Although the BR curve "looks" like a straight line 
on a log-log graph (See Fig. 3) for smallest, we cannot entirely 
reconcile our result with those of Larsen (1977) . However, 
we note that the suggested BR curve reduces to the form given 
by O'Gara (Heck and Brandt, 1977) for plant response. This 
can be demonstrated as follows: For small o(t (3 J can be 
written as 



r n (9b) 

^ r Ut 2~ 

Equation (9b) is of the form given by O'Gara who obtained it by 
fitting a curve to his experimental data. 

We realize that the determination of al and C for 
various pollutant-receptor pairs presents major difficulties. 
The importance of these parameters is tied to the validity of 
our model. However, even if the model is not realistic the 
principles behind the derivation of air quality criteria and 
our air pollution index have importance of their own. 
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